Doping incorporation paths in catalyst-free Be-doped GaAs nanowires 
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The incorporation paths of Be in GaAs nanowires grown by the Ga-assisted method in molecular beam epitaxy 
has been investigated by electrical measurements of nanowires with different doping profiles. We find that 
Be atoms incorporate preferentially via the nanowire side facets, while the incorporation path through the 
Ga droplet is negligible. We also demonstrate that Be can diffuse into the volume of the nanowire giving an 
alternative incorporation path. This work is an important step towards controlled doping of nanowires and 
will serve as a help for designing future devices based on nanowires. 



Semiconductor nanowires have stimulated extensive in- 
terest in the last decade because of their potential as 
building blocks in future generations of electronic, op- 
toelectronic devices as well as for energy conversion and 
applications 1 5 . In order for nanowires to become a tech- 
nological reality, control of the conductivity by doping 
is extremely important. Doping involves the incorpora- 
tion of impurities with a small ionization energy, which 
can transfer carriers to either the conduction or valence 
band. The incorporation paths of dopants in nanowires 
have been discussed in the last few years, in views of con- 
trolling their concentration and position in the nanowire. 
For this, the understanding of the growth process of 
a nanowire is important. Bottom- up grown nanowires 
are typically formed via the Vapor-Liquid- Solid (VLS) 
mechanism 6 8 , where a nanoscale liquid droplet acts as 
a catalyst for nanowire crystal formation. It is generally 
accepted that incorporation paths of dopants are sim- 
ilar from the growth precursors and result either from 
the radial growth of the nanowire or through diffusion 
through the catalyst 9 14 . Still, one may consider a third 
path, which consists on the diffusion of dopants in the 
core and/or from the shell to the nanowire core. The 
relevance of this path should depend on the coefficient 
of diffusion of the dopant and the growth temperatures 
used 15 . Typical dopants with a high diffusion coefficient 
in III-Vs are Be and Te 16,17 . In this work, we investigate 
the doping of GaAs nanowires with Be. By measuring 
the spatial dependence of the conductivity of nanowires 
deposited with a flow of Be under different conditions we 
identify the incorporation paths and discuss the advan- 
tages and limitations of this dopant. In principle, this 
work can be extrapolated to other III-V materials. 

In the correlation between the electrical conductiv- 
ity and dopant incorporation, three main incorporation 
paths will be distinguished as schematized in Fig 1): i) 
Axial incorporation through the VLS mechanism, ii) ra- 
dial incorporation through the VS mechanism and as an 



extension of the latter and ii) diffusion of the dopants 
from the shell to the core in the volume of the nanowire. 
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FIG. 1. Processes that influence the Be incorporation in GaAs 
nanowires: vapor- liquid-solid mechanism (VLS), possibility of 
growing a doped shell (VS) and diffusion of dopants during 
the growth process. 

The nanowires were grown on (111) Si using a self- 
catalyzed (VLS) method 18 20 , with a Ga deposition rate 
corresponding to a nominal growth of 0.27 A/s , for times 
ranging between 30 and 60 minutes (Tab I), at 630° C 
substrate temperature and a V/III flux ratio of 60 21 . 
The vertical growth rate for all the nanowires was around 
15 fim/h. The p-doping was achieved by adding a flux of 
beryllium either during growth of the core or a posteriori 
during the growth of a shell. The shell was obtained at a 
lower temperature (465° C), switching the As source from 
AS4 to As2 and increasing the V/III ratio to 150 22 . To 
give some insight into the Be incorporation mechanisms 
during axial VLS and radial VS growth of self-catalyzed 
GaAs nanowires, six different types of growths with vary- 
ing doping profiles were carried out, see Table I. 

To give a reasonable estimate of the effective carrier 
concentrations, the effective resistance in the nanowires 
was measured by carrying out 4-point electrical measure- 
ments at room temperature 23 on around 50 samples. The 
high number of devices were obtained within a reasonable 
time frame with our auto-contacting software 24 . The 
electrical contacts consisted of Pd/Ti/Au (40/10/250 
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5 • 10 18 


0.06 ± 0.03 
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60 





2.5 • 10 19 


16000 ± 5000 
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45 


30 


1 • 10 18 


3±2 
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45 


30 


5 • 10 18 


50 ±20 


5 


30 


10 


1.5 -10 19 


1000 ± 600 


6 


30 


30 


1.5 -10 19 


2200 ± 500 



TABLE I. Time of axial nanowires growth, nominal shell 
thickness, nominal doping concentration (corresponding to 
the planar growth under the same conditions) and measured 
conductivity obtained by 4-point contact configuration. 



nm). In order to ensure ohmic and reproducible contacts, 
the contacted samples were annealed at 300°C for 20 min- 
utes. Figure 2a) shows a typical contacted nanowire for 
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FIG. 2. a) I-V examples of each growth performed with sin- 
gle 2 contact configuration. The linearity of the curve shows 
that Pd/Ti/Au electrical contacts are ohmic on Be doped 
GaAs nanowires. In the insert a SEM image of a contacted 
nanowire. b) Nanowires conductivity obtained from 4 point 
measurements and a corresponding SEM image. 

single 2 point contact configuration. Multiple contacts 
were also realized for understanding the spatial depen- 
dence of the conductivity along the nanowire axis. Exam- 
ples of 2-point current- voltage characteristics are shown 
in the main graph of Fig. 2a). The linearity of the curves 
indicates that all contacts are ohmic, even for the highly 
resistive samples. The conductivity of the nanowires ob- 
tained from 4-point probe measurements can be found in 
Tab I and in Fig. 2b). The apparent conductivity ob- 
tained in the nanowires grown with a nominally doped 
core and without a shell (sample 1) is much lower as 
compared to the nominal doping concentration. Sample 
4 has been grown with the same core conditions as Sam- 
ple 1 with an additional shell of 30 nm with the same 
Be flux. By comparing the two samples we can sepa- 
rate the transport contribution due to the VLS step and 
obtain the contribution to the conductivity from the VS 
step. The doping contribution from the VLS step is neg- 
ligible indicating that the flux of Be atoms to liquid-solid 
growth front is relatively small 25,26 . We therefore deduce 



that the incorporation path of Be through the droplet can 
be neglected. 

The average carrier concentration in GaAs at room 
temperature can be obtained from the relation p = Na, 
where Na is the doping density. One can use the equa- 
tion: 



N A 



(i) 



where e is the electrical charge of the electrons, \i is 
the mobility in the nanowires and a is the conductivity. 
We assume /i = 31 cm 2 /(V • s), as recently measured 
in similar nanowires and doping concentration range 27 . 
The fundamental parameter we extract from the electri- 
cal measurements is the current crossing the nanowire. In 
order to extract the doping concentration, the following 
effects should be taken into account: (i) non-uniform ra- 
dial doping distribution, (ii) surface effects such as deple- 
tion and dopant deactivation and (iii) diffusion of dopants 
during the growth process. Our model considers these 
three effects, as it is shown in the following. 

We start by presenting the effect of surface band bend- 
ing, especially important for small diameters and/or low 
doping concentrations. The exposure of nanowires to am- 
bient conditions causes the formation of a thin oxide layer 
on the surface. This results in the pinning of the Fermi 
level at the surface. In GaAs, the pinning occurs near 
the middle of the bandgap, thereby producing a depletion 
layer close to the surface 28 . In a semi-classical model the 
spatial dependence of the band bending is obtained by 
considering the difference cp of the Fermi level at surface 
states with respect to the Fermi level of the bulk. The 
variation of cp in the nanowire geometry can be described 
by the Poisson equation in cylindrical coordinates: 



1 d dip(r) 
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This equation can be solved by using the boundary 
conditions of a vanishing electric field at the surface: 
^(^o) = $7 where corresponds to the pinning posi- 
tion of the Fermi level (0.5V in p-type GaAs 29 ). The 
result is an implicit equation for the nanowire depletion 
layer width w which depends on the nanowire radius ro 
and the doping concentration Ngheii 30 '- 
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The values of the depletion region w as a function of the 
doping concentration and for two different nanowire radii 
is shown in Figure 3. For high doping concentrations, the 
depletion region can be as small as few nm. It increases 
rapidly for lower doping concentrations. Already in the 
case of doping concentration close to 10 18 atoms / cm 3 , the 
depleted region corresponds to several tens nm. For ex- 
ample, a nanowire with ro = 40 nm would be completely 
depleted and therefore insulating for carrier concentra- 
tions below 1 • 10 18 atoms / cm 3 . 
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FIG. 3. Calculated depletion region width in dependence of 
nanowire concentration for radii ro = 40 nm and ro = 100 
nm. 



The actual doping profile in the nanowire depends also 
on the diffusion process during growth, which is non- 
negligible for Be in GaAs 15 ' 31 . In fact, during axial 
growth there will be a concentration of Be adatoms on the 
nanowire sidewalls which can be incorporated thanks to 
the non-zero radial growth rate. In fact, according to pre- 
vious studies on similar nanowires 18 ' 23 , the radial growth 
rate is about 1000 times lower than the axial growth rate. 
The concentration of Be in the vicinity of the nanowire 
surface is kept constant due to the constant incident Be 
flux, which leads to a quasi steady state concentration at 
the surface during growth, po. Diffusion of Be from the 
surface to the nanowire core is driven by the gradient in 
Be concentration. The diffusion is given by 32 



p(x, t) = po • erfc 



2-(Dt) 1 / 2 



(4) 



where po is the doping concentration at the interface, x 
is the distance from the interface, t is the diffusion time 
and D is the diffusion coefficient 



D = D -e 



-E /kT 



(5) 



where for Be diffusion in GaAs, Z^o = 0.655cm 2 /s and 
E = 2.43eV 32 . The diffusion length is defined as the 
distance with which the concentration is 1/e of the shell 
concentration (Table II). As schematically drawn in Fig. 
4, the existence of a depletion region and the diffusion of 
Be strongly modify the range of electrically active part 
of the nanowires. 

Since the 4-point measurements are performed in the 
center of the nanowire, it is assumed that the diffusion 
at that point corresponds to half of the full VLS growth 
time. The full nanowire cross-section is integrated and 
the depletion region included as described by equation 
(3). Then, the doping concentration and electrical carrier 
concentration can be deduced from the electrical mea- 
surements: 
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FIG. 4. a) Section view of a nanowire. The doping and the 
carrier concentration are reported as a function of position 
along the nanowire diameter, b) Comparison between carrier 
density concentration (red) and calculated doping concentra- 
tion in the shell (black). The blue dotted line represents 
N a = N s heii- The error bar reported for every dots repre- 
sent only the standard deviation calculated on the electrical 
measurements. 



where R/L is the resistance per unit length, N s h e ii is the 
shell concentration and the doping distribution (p(x,t)) 
is: 



p(x,t) 
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(1 ±0.5) • 10 15 


(1.2 ±0.6) • 10 18 
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38 


2 


(3± 1) • 10 19 


(5±1) • 10 19 


8 


3.5 


53 


3 


(5±3) • 10 15 


(5±1)-10 17 


36 


50 


46 


4 


(1 db 0.5) • 10 17 


(6 ±3) • 10 17 


36 


39 


46 


5 


(3 ± 2) • 10 18 


(5 ±3) • 10 18 


14 


12 


38 


6 


(4 ± 1) • 10 18 


(7 ±2) • 10 18 


34 


10 


38 



TABLE II. Carrier concentration, shell doping concentration 
calculated with the model described previously, total expected 
shell, depletion region (w) and diffusion length defined by the 
distance which the concentration is 1/e of the shell concen- 
tration. 

For doping concentrations below 5 • 10 18 cm -3 the de- 
pletion region is larger than the shell thickness. This 
concerns sample 1,3 and 4. These three samples exhibit 
a very low conductivity. By taking into account the de- 
pletion region at the surface, we can extract the dop- 
ing concentration at the shell and find that it is very 
close to the nominal doping concentration. For higher 
doping concentrations, the depletion region reaches few 



nm and affects much less the overall conductivity. We 
should point out that electrical measurements performed 
along the nanowire axis with multiple contacts showed 
a homogeneous conductivity along the nanowire. It is 
interesting to note that even though the incorporation 
of Be via VLS mechanism is negligible, it is possible to 
dope the nanowires almost homogeneously during axial 
growth provided the Be flux is high. For nominal doping 
concentration higher than 10 19 cm -3 the highly doped 
shell and the long growth time results into a strong dif- 
fusion of Be inside the core, leading to an almost fully 
doped nanowire 33 . This is extremely advantageous for 
example in the case where an ohmic contact with a sub- 
strate is relevant (e.g. solar cells). 

In conclusion, we have shown that the Be atoms are 
mostly incorporated from the side facets and that the 
incorporation through the Ga droplet is negligible. The 
doping concentration is homogeneous along the nanowire 
and can be tuned between 6 • 10 1T and 5 • 10 19 cm -3 . 
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